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Dwarf spheroidal galaxies have a large mass to light ratio and low astrophysical background, and 
are therefore considered one of the most promising targets for dark matter searches in the gamma-ray 
band. By applying a joint likelihood analysis, the power of resultant limits in case of no detection 
can be enhanced and robust constraints on the dark matter parameter space can be obtained. We 
present results from a combined analysis of 10 dwarf spheroidal galaxies using Fermi-LAT data. 
Different annihilation channels have been analyzed and uncertainties from astrophysical properties 
have been taken into account. 



I. INTRODUCTION 

The Fermi Gamma-ray Space Telescope was 
launched on June 11, 2008. Its main instrument, the 
Large Area Telescope (Fermi-LAT), observes the en- 
tire sky every ^3 hours (2 orbits) with a field of view 
covering ~2.4 sr and a sensitive energy range extend- 
ing from 20MeV to >300GeV PQ. These properties 
make the Fermi-LAT an excellent instrument for dark 
matter (DM) searches. 

One of the leading DM candidates is a weakly in- 
teracting massive particle (WIMP). The gamma-ray 
flux from self-annihilating WIMPs can be expressed 
as 4>wimp(E,iP) = J(ip) x <& PP (E) (see e.g., [2], see 
also [3] for a review), where <& PP (E) is the "particle 
physics factor" described by 
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and J(V0 is the "astrophysical factor", or J-factor, 
described by 
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Here (av) is the (velocity-averaged) WIMP annihi- 
lation cross section times relative velocity, mwiMP is 



the WIMP mass, Ny/(E) is the gamma-ray energy dis- 
tribution per annihilation, and p(r) is the dark matter 
density distribution. 

Dwarf spheroidal galaxies (dSphs) are considered to 
be DM dominated systems since they have a very high 
mass to light ratio. They are near-by (~ lOOkpc) and 
they have low background since most dSphs are ex- 
pected to be free from other astrophysical gamma-ray 
sources and they have a small gas content. This makes 
them interesting targets for gamma-ray DM searches. 
However, the expected flux from DM annihilation or 
decay for dSphs are expected to be very low. The 
Fermi-LAT collaboration has recently presented re- 
sults from a DM search in a number of dSphs [I]. 
In this work we perform a joint likelihood analysis, 
taking advantage of the fact that the DM spectra are 
the same in all targets. 

We present results considering ten dSphs taking 
into account the uncertainties in the astrophysical fac- 
tors. The results presented in these proceedings are 
updated with respect to the preliminary results pre- 
sented at the conference. This work was initially pre- 
sented in [5], and final results are presented in [6]. 



II. ANALYSIS 
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We have observed ten dSphs, listed in Table|TJ which 
is the same set of dwarfs for which annihilation cross 
section limits were presented in [4] with the addition 
of Carina and Segue 1, using 24 months of data. We 
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TABLE I. Position, distance, and J-factor (under assump- 
tion of a Navarro-Frenk- White profile) of each dSph. The 
4th column shows the mode of the posterior distribution of 
log 10 J, and the 5th column indicates its 68% C.L. error. 
See the text for further details. The J-factors correspond 
to the pair annihilation flux coming from a cone of solid 
angle AQ — 2.4 • 10 -4 sr. The final column indicates the 
reference for the kinematic dataset used. 



Name 
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d lo 


gio(J r ) 


a 


ref. 




deg. 


deg. 


kpc lo 


gioIGeV' 


! cm" 5 ] 




Bootes I 


358.08 


69.62 


60 


17.7 


0.34 


PI 


Carina 


260.11 


-22.22 


101 


18.0 


0.13 


[7] 


Coma Berenices 


241.9 


83.6 


44 


19.0 


0.37 


9 


Draco 


86.37 


34.72 


80 


18.8 


0.13 


.7 


Fornax 


237.1 


-65.7 


138 


17.7 


0.23 


7 


Sculptor 


287.15 


-83.16 


80 


18.4 


0.13 


7 


Segue 1 


220.48 


50.42 


23 


19.6 


0.53 


8. 


Sextans 


243.4 


42.2 


86 


17.8 


0.23 


7 


Ursa Major II 


152.46 


37.44 


32 


19.6 


0.40 


9. 


Ursa Minor 


104.95 


44.80 


66 


18.5 


0.18 


7 



have used the diffuse event class which only contains 
the events with the highest gamma-like confidence, 
and we have chosen events ranging from 200 MeV 
to 100 GeV. We used the Fermi-LAT instrument re- 
sponce function P6_V3 -DIFFUSE. Our region of inter- 
est (ROI) is a region of 10 degrees radius centered on 
dSph location. Standard cuts removing Earth albedo 
photons have been made. The dSphs are modeled as 
DM point sources using the DMFit package [11] where 
we consider 100% annihilation into the bb, the t + t~, 
the W + W~ , and the fi + [i~ annihilation channel. The 
background is modeled according to Fermi-LAT rec- 
ommendations |12j . and sources within 15 degrees are 
modeled according to the first year point source cat- 
alog [15] . We perform a binned analysis to use both 
energy and spatial information. The data selection 
and analysis are performed using the Fermi-LAT anal- 
ysis package, ScienceTools [13] . and the upper limits 
are obtained using profile likelihood as implemented 
in the MINUIT processor MINOS [H]. 

One large uncertainty in indirect DM detection 
methods arises from the uncertainties in the astro- 
physical factors, the J-factors. We have included these 
uncertainties by including the distribution of the J- 
factors in the likelihood fit, treating them as nuissance 
parameters. This is the first time J-factor uncertain- 
ties are included in this way. In our fits, the parame- 
ter of interest is the WIMP annihillation cross-section, 
{av), and the nuissance parameters are the J-factors, 
the normalizations of the Diffuse Backgrounds and 
the normalizations of sources within 5 degrees of the 
dSphs. 

With this addition, the joint likelihood considered 



in our analysis becomes: 
L(D|pw,{p}i) =n^ AT ( £) IPw,Pi) 

i 

x 1 c -(log 1n (./ t )-log 1f ,(J I )) 2 /2 CT ? 

ln(10) J l V^a l 

(3) 

where L\ AT denotes the binned Poisson likelihood 
that is commonly used in a standard single ROI analy- 
sis of the LAT data, i indexes the ROIs, D represents 
the binned gamma-ray data, pw represents the set 
of ROI- independent DM parameters ((av) and mw), 
{p}i are the ROI-dependent model parameters. In 
this analysis, {p}i includes the normalizations of the 
nearby p oint and diffuse sources and the J-factor, J, ; . 
log 10 (Ji) and cr,j are the mean and standard deviation 
of the distribution of log 10 (J,), approximated to be 
gaussian. The values for the J-factors have been up- 
dated since the preliminary work was presented at the 
conference, and the updated values are found in Table 
[i] A detailed discussion about how the J-factors are 
obtained are found in [6]. 

We have performed tests on the coverage of the joint 
likelihood method and also tests on how the combined 
limits scale with the number of added targets. In our 
fits, we constrain the (av) parameter to be positive. 
Our tests show that this will make the combined limit 
scale as better than l/v^V, where N is the number 
of added targets. The hardness of the DM spectra 
will also influence, resulting in a larger improvement 
for harder spectra. The coverage of our method is 
good but with slight overcoverage for small signals, 
i.e. conservative limits. 



III. RESULTS 

In Fig. [T]we present constraints on the WIMP anni- 
hilation cross-section into 100% bb for ten dSphs, both 
individual and combined limits. In Fig. [2] we present 
the combined limits on WIMP annihilation cross sec- 
tion for annihilation into 100% bb, 100% t+t', 100% 
W + W~ , and 100% • J-factor uncertainties have 

been included. The joint likelihood method allows us 
to rule out WIMP annihilation with cross sections pre- 
dicted by the most generic cosmological calculation up 
to a mass of ^30 GeV for the bb and the t + t~ anni- 
hilation channel. 

Final results are presented in [6]. 
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FIG. 1. 95% Upper limits on WIMP annihilation cross sec- 
tion for annihilation into 100% bb. The expected thermal 
WIMP cross-section is plotted as a reference. 
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FIG. 2. 95% Upper limits on WIMP annihilation cross 
section for annihilation into 100% bb, 100% t+t~, 100% 
W + W~, and 100% The expected thermal WIMP 

cross-section is plotted as a reference. 
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